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We study the diffusiophoretic self-propulsion of a colloidal catalytic particle due to a surface 
chemical reaction in a vicinity of a solid wall. Diffusiophoresis is a chemico-mechanical transduction 
mechanism in which a concentration gradient of an interacting solute produces an unbalanced force 
on a colloidal particle and drives it along the gradient. We consider a spherical particle with 
an axisymmetric reacting cap covering the polar angle range 0 < ^ < Ocap in the presence of a 
repulsive solute, near an infinite planar wall, and solve the coupled solute concentration and Stokes 
equations, using a mixture of numerical and analytic arguments. The resulting particle trajectory 
is determined by Ocap and the initial orientation of the symmetry axis with respect to the plane. At 
normal incidence the particle initially moves away from or towards the wall, depending on whether 
the cap faces towards or away, respectively, but even in the latter case the particle never reaches the 
wall due to hydrodynamic lubrication resistance. For other initial orientations, when Ocap < 115° 
the particle either moves away immediately or else rotates along its trajectory so as to cause the 
active side to face the wall and the particle to rebound. For higher coverage we find trajectories 
where the particle skims along the wall at constant separation or else comes to rest. We provide a 
phase diagram giving the nature of the trajectory (repulsion, rebound, skimming or stationary) as 
a function of Ocap and the initial orientation. 


* These two authors contributed equally to this work. 
^ cmaldarelli@ccny.cuny.edu. 



2 


I. INTRODUCTION 


The newest generation of ultra-miniaturized engines now under development are constructed using colloidal sized 
objects, which are engineered to react with solute in a solution in which the colloids are immersed, and convert the 
chemical reaction energy into mechanical self-propulsion without moving parts (for reviews see Howse [1], Sen and 
Velegol [2-4], Wang [5-9], Schmidt [10], Pumera [11, 12], Kapral [13], and Ozin [14]). These autonomous molecular 
locomotors or artificial swimmers represent a bottom-up approach to mechanical design made possible by advances 
in micro and nano-fabrication, and are at the center of a wide range of potential applications at the micro- and nano¬ 
scale. The locomoters themselves can be used as micropumps [15], or as roving sensors if their chemical propulsion 
is affected by the presence of solute [16]. The locomoters can be configured with chemical receptors to function as 
transporters for capturing, towing and delivering molecular cargo [9, 17, 18] as for example drugs to selected targets 
[19-21], molecular building blocks to assemble supramolecular structures [22], or analytes in microfluidic networks for 
lab on a chip operations [23-26]. Colloidal motors with chemical receptors have also shown promise as the engines of 
shuttles for the capture and movement of biological cells [27-29] . 

The operating principle for the self-propulsion of molecular motors is the chemico-mechanical transduction mech¬ 
anism which converts the reaction energy into the propulsive motion. The first set of synthetic locomotors were 
bimetallic rods made of platinum and gold (a few microns in length and a few hundred nanometers in diameter) pre¬ 
pared by sequential electrodeposition in a membrane template. A catalytic oxidation of a solute, hydrogen peroxide 
takes place on the platinum (anode) side of the rod to produce protons in solution, and electrons. The protons migrate 
to the opposite gold (cathode) side of the rod through the solution, while the electrons conduct through the rod to the 
gold side. At this cathode end, hydrogen peroxide combines with the released protons and the electrons conducted 
through the metal to produce water [30-36]. The accumulation of protons at one end of the bisegmented rod, relative 
to the opposite end creates an electric dipole field around the rod directed from the anode to the cathode. The field 
acts on the diffuse layer of positive charge in the water around the cylindrical sides of the particle (which balances 
the negative charge of the metal) driving fluid from the anode to the cathode and propelling the rod in the opposite 
direction. 

More recent research has adopted a tubular configuration for the engine, rather than a nano-rod structure, in order 
to increase the efficiency of the hydrodynamic propulsion. Using either thin-film roll-up technology [10, 37-39], or 
template membrane electrodeposition, conically shaped hollow cylinders are formed tens of microns in length and a 
few microns in diameter with an inner catalytic platinum surface which catalyzes the reaction of hydrogen peroxide 
to water and oxygen on the inside of the cylinder. The buildup of oxygen creates bubbles which are expelled through 
the wider opening of the cylinder, creating a recoil which propels the object. 

A third molecular engine configuration, which we will focus on in this study, is constructed around spherical 
colloids, and has the advantage that it is the easiest of the engines to fabricate. In this design, as shown in Fig. 1, 
nonconducting (usually polystyrene) spheres a few microns in diameter are first coated on one side with platinum to 
form a Janus colloid (see Howse, Golestanian and coworkers and references [1, 40-42]). The functionalized colloids 
are then immersed in a hydrogen peroxide solution; because the underlying polystyrene is nonconducting, there is no 
electrochemical driving force for oxidation, and the platinum catalyzes the hydrogen peroxide in solution directly to 
water and oxygen. The experiments are usually undertaken in two regimes: The first is at a high hydrogen peroxide 
concentration, in which oxygen production exceeds the equilibrium solubility in aqueous solution, and microbubble 
production is observed on the active side of the colloid. As in the tubular jet configuration, the recoil force on 
the colloid projects the particle in the direction of the passive side. The second case is one in which the hydrogen 
peroxide concentration is low, and oxygen bubbles are not evident, but the colloid is still observed to move in the same 
direction (leading with the passive side). Howse, Golestanian and collaborators [40, 41] have argued that unbalanced 
intermolecular interactions between the hydrogen peroxide reactant and the oxygen product create the propulsion. 
This type of motion can be explained by the diagram in Fig. 1(a). Assuming that the (oxygen) product species 
produced by the reaction repels the colloid particle, the excess of this product on one side of the colloid due to the 
asymmetric distribution driven by the catalytic conversion leads to a net repulsive force projecting the particle in 
the direction of the passive side. The intermolecular interaction occurs over a small length scale, denoted by L in 
the figure, and is of the order of a few nanometers. This propulsion force is termed diffusiophoresis after Derjaguin 
[47, 48] and later Anderson et al [43-46] who studied this problem for the case of a gradient in solute applied directly 
across the particle (rectified diffusiophoresis). Although, the neutral solute gradient driven self-diffusiophoresis has 
been widely considered as an accepted mechanism for the self-propulsion of Janus catalytic swimmers due to hydrogen 
peroxide decomposition on the platinum surface, recent experimental studies reveal that the electrochemical effect 
due to ion diffusion can not be entirely neglected and further investigation should be undertaken to elucidate and 
disentangle these effects [49]. 

Unbalanced intermolecular interactions which are caused by an asymmetric reaction around a particle, and which 
in turn drive a propulsive motion, represents a general type of mechanism by which a molecular engine can operate. 
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FIG. 1. A Janus colloid molecular motor: (a) A difference in concentration of product (oxygen) created by the reaction of hydrogen 
peroxide on the catalyst half of a Janus particle generates an unbalanced intermolecular force on the particle (shown here as repulsive 
and acting over a length scale L) which propels the particle towards the lower concentration, shown by the particle velocity vector V. (b) 
The thin film hydrodynamic analysis of Anderson et al [43-46] for a ^ L: Greater net repulsion of the solute against the particle on the 
high concentration side in the intermolecular interaction layer generates a lower pressure Pi relative to outside the intermolecular layer 
(Pouter) to balance this repulsion compared to the pressure P2 at the low concentration side. This causes a streaming fiow towards the 
high concentration side (a slip velocity) forcing the particle to moves in the opposite direction or leading with the passive side as in (a). 


and is adopted in this study. For particles in an infinite liquid phase hydrodynamic models have been developed to 
predict the diffusiophoretic velocity [41, 50-53]. The hydrodynamics was first studied theoretically in a general way 
by Derjaguin [47, 48] and later Anderson et al [43-46] in the case in which gradients of solute are applied directly 
across the colloid. These authors studied the problem in the limit in which the interaction length L is much smaller 
than the particle radius a, i.e. A = L/a <C 1, and the particle moves through a continuous Newtonian phase in Stokes 
flow. Within a continuum hydrodynamic framework, the effect of the intermolecular interaction of the particle with 
the solute is accounted for by adding a body force to the Stokes equations, —CV0, where C is the solute concentration 
and (j){r) is the potential energy of interaction of the solute with the colloid [43, 54]. In the limit A = L/a ^ 1, the 
flow field divides into a thin inner region (within a distance L of the particle much smaller than the radius a ) where 
the interaction force is effective, and an outer region where it is relatively unimportant (Fig. 1(b)). The inward 
repulsion of the solute with the particle creates a lower pressure in the thin layer (relative to outside the layer) which 
is necessary to balance this outward force. Since the repulsion of the solute to the particle is greater along the side of 
the particle where the solute concentration is larger, the pressure is lowest there, and this generates a streaming flow 
to the higher concentration side as the particle moves in the direction of lower concentration (towards the passive 
side). This streaming flow creates a tangential velocity VsUp at the outer edge of the inner region which is a function 
of the concentration gradient along the surface at this edge, VsUp = —hVsCs where s measures distance along the 
surface, Cg is the concentration of solute along the outer edge and 6 is a slip coefficient which is a function of the 
intermolecular forces and the concentration distribution of solute in the inner region. From the perspective of the 
outer region, the flow at the outer edge of the inner region represents a “slip” velocity on the particle surface. To 
leading order, the flow in the outer region and the translational diffusiophoretic velocity V of the particle are solved 
in terms of this slip velocity. The solution of the mass conservation equation for the concentration field in the outer 
region allows for the calculation of the gradient in the solute concentration at the outer edge of the inner region from 
which the particle velocity obtains. Anderson et al [43-46] found for the slip coefficient 
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b = -/ y [exp {-4>{y)/kBT) - 1] dy, (1) 
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where 0(r) is the interaction potential energy, ksT is the thermal energy, /i is the solution viscosity and y = {R—l)/X 
is the rescaled radius with A = L/a. More recently, the validation of this continuum framework for sub-micron size 
colloidal particle has been tested against molecular dynamics simulation study [55] and a micro-mechanical colloidal 
perspective [56]. 

The self-diffusiophoresis of a non-Brownian Janus colloid with thin interaction layer and a well-defined net potential 
interaction between solute and the colloid in an infinite medium has been now well-established and it was shown that 
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the swimming velocity of the colloid can be described via two non-dimensional numbers so called Damkohler Da, and 
Peclet Pe numbers (see below for their definitions). In order to describe the hydrodynamics of self-diffusiophoresis, 
in which the gradient responsible for driving the phoretic motion is generated by a surface reaction, the solutal 
concentration fields have to be resolved for the reaction-diffusion problem. The hydrodynamic studies assume an 
irreversible surface reaction of a reactant to products which occurs only along the reactive cap (Fig. 1(b), 0 < 6> < Ocap)^ 
and produces a constant flux (dimensional) Nq of product. Along the rest of the colloid surface, the solute fluxes 
are zero since the solutes can not penetrate into the particle. This assumption of constant flux along the reactive 
part is equivalent to assuming the reaction rate is slow relative to diffusion of the reactants [54] and product (small 
Damkohler {Da) number where Da is defined generally as Da = ^ where k is the surface rate of reaction of the 
reactant A and D in the diffusion coefficient of the reactant. In the limit Da ^ 0, to leading order, the concentration 
of the solutes is spherically symmetric, and to first order, the flux off of the cap, due to the reaction, is a constant 
and of order Da. To calculate the self-diffusiophoretic velocity V, the hydrodynamic studies assume Anderson’s 
expression for the slip velocity (Eq. 1). To compute the concentration field at the outer edge of the inner region, 
two approximations are made. First, only the concentration field in the outer region is solved, with the constants 
in the outer solution evaluated by applying the boundary conditions directly at = 1 (the colloid surface). This 
computation then provides the solution for 1,0) = Cs{0). As a second approximation, the Anderson 

expression for the slip velocity coefficient (Eq. 1) is used to calculate the velocity. With these assumptions, the studies 
indicate that for the case in which only the product “B” interacts with the colloid, the self-diffusiophoretic velocity is 
to leading order in A given by: 


T/ 

Voo-h D 


cos^6>, 


cap 


( 2 ) 


where Wo is the dimensional diffusiophoretic velocity. For a fixed value of the dimensional flux, Nq, the velocity 
is again found to be independent of the particle radius a. The above expression has been used by Golestanian et 
a/[41, 50, 51] to model the observed motion of Janus colloids driven by the asymmetric decomposition of hydrogen 
peroxide to oxygen on the platinum functionalized side of the particles. Note that the above relationship assume 
the colloid motion is non-Brownian, i.e. the time scale for rotational diffusion of the colloid is much longer than 
the time scale for the directed motion of the particle ^src^nian ^ 2 , the typical value of this ratio for micron-size 

Uc 

Janus colloid which swim with the speed of 5 — 50^ is about 30 — 6000 and therefore it is a legitimate assumption 
to ignore Brownian rotation of the colloid in the calculation for the short time t ^Brownian, however the role of 
Brownian rotational diffusion will be important in long time and consequently results in a diffusive behavior [50]. 
Later on, using the same Anderson mobility framework, Popescu et a/. [52, 53] extended the above expression for 
colloid particles with axisymmetric, ellipsoidal shapes. For intermediate to high values of Da, it was shown that the 
swimming velocity of the colloid increases monotonically until it reaches a plateau where the system is diffusion limited 
[57]. The influence of solute advection for rectified diffusiophoresis was studied by Keh and Weng [58] and recently 
Khair [59]. It was found that for a thin interaction layer compared to colloid size, the effect of solute advection to 
the diffusiophoretic velocity is O(Pe^) and particle velocity persistently decreases as Peclet number increases (Peclet 
number is the ratio of solute advection to solute diffusion which defined as ^ where Uc is the characteristic 
diffusiophoretic velocity). In self-diffusiophoresis on the other hand, Milchen and Lauga [60] found that the effect of 



FIG. 2. A Janus swimmer with active area of size 6cap near a solid boundary, with active area (in red) tilted at an angle S from the 
perpendicular to the boundary surface and with the colloid a distance S from the wall. S = 0 represents a colloid with its active area 
directly facing the wall. 
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solute advection on diffusiophoretic velocity is 0{Pe) and hence depending on attractive or repulsive net interaction of 
the solute with the colloid, the solute advection can increase or decrease the self-diffusiophoretic velocity respectively. 
More recently, the effect of external shear flow in self-diffusiophoresis has been examined by Frankel and Khair [61]. 
The thin interaction layer assumption and the slip velocity argument can be relaxed by introducing a body force 
term in Stokes equation which is a consequence of net interaction between the solutes and also a new term would 
appear in mass balance conservation which is a contribution to the flux of solute due to the net interactions with the 
colloid. In this case, it was shown that the swimming velocity of the colloid decreases as the size of interaction layer 
becomes larger [54, 62]. An alternative approach to study the colloidal particle motion in the presence of gradient 
of another species is the micro-mechanical perspective in which the particle and solutes are Brownian particles in 
continuum solvent and in this treatment the constraint of small and low volume fraction of solutes can be relaxed 
[56]. Moreover it was shown the micro-mechanical perspective can successfully recapitulate the continuum treatment 
of diffusiophoresis in the region of overlap [54, 56]. The shape of the particle itself can also brings about propulsion 
of the particle, in this case it was shown that a non-spherical particle is able to move even with uniform catalytic 
activity [63, 64]. What remains still as a potential challenge is a controlled motion at micro-nanoscale and it’s been 
suggested that by applying an external field, the Janus swimmers can be steered towards a specific location [28]. 

The investigations discussed above are an accounting of the self-propulsion of a reactive Janus swimmer in an 
infinite medium. Since a lot of anticipated applications of colloidal motors involve in the motion of these swimmers in 
the vicinity of an interface or a solid wall, it’s crucial to understand the effect of interface on the dynamics of colloidal 
particle. Moreover, we want to respond the question weather it is feasible to steer these artificial microswimmers 
via an interface similar to the observation for biological swimmers, e.g. spermatozoid [65] or bacteria [66] . Despite 
swimming in the vicinity of a boundary can display a strikingly more complex set of trajectories than the rectilinear 
motion the swimmer executes in an infinite fluid space. The richer dynamics of colloidal propulsion near a boundary 
derives from the diffusive interaction of the solute concentration field generated at the surface of the active area of the 
colloid with the boundary which changes the diffusiophoretic propulsive force. In addition, there is related work on 
electrophoretic and dipolophoretic motion of colloids [67-70] and biological swimmers such as bacteria [71, 72] near 
walls and finally the motion of biological swimmers, e.g. bacteria close to interfaces is crucial for understanding in 
broad range of problems involving biofilm formation on surfaces [73-76]. 

In this study we will consider an infinite solid surface as the boundary, which we assume cannot be penetrated by the 

solute (see Fig. 2). For a colloidal particle which is self-diffusiophoretically swims in a proximity of a solid wall, the 

role of the wall is not only important because of the hydrodynamic interactions, but also it distorts the concentration 

field around the particle and consequently alters the dynamics of colloidal motors. We denote by S the tilt angle of 

the outward normal of the center of the spherical cap of the active area of the swimmer with the perpendicular to the 

boundary (defining S = 0 to be the orientation in which the active area normal points towards the boundary), and 

S the edge-to-edge perpendicular distance of the sphere to the wall. We will assume Stokes flow, and a slow reaction 

so that the flux of solute on the active area is specified. The mass transfer is regarded as quasi-static, i.e. the time 

2 

scale for diffusive equilibration ^ is much shorter than the time scale for movement of the particle over its radius a, 
with the velocity of the particle as estimated by the velocity far from the boundary, Voo^ Eq. 2. This corresponds 
to Peclet numbers Pe = 1^ where D is solute diffusion coefficient. In addition, we will also assume that the 

intermolecular interaction length scale L is much smaller than the colloid radius a so that the simplifying framework 
in which propulsion is driven by a slip velocity on the surface dictated by the gradient in the surface concentration of 
the solute around the colloid. As in our discussion above, we will assume that the catalytic reaction on the surface 
produces a solute which acts repulsively with the colloid so that the accumulation of solute on the active side of 
the swimmer pushes the swimmer in the opposite direction, or equivalently creates a slip on the active area in the 
direction away from the passive side. Hence the propulsion velocity V (relative to the value far from the wall, Eoo? 
eq. 2) becomes an instantaneous function of the wall-colloid separation distance, J, the tilt S and the size of the 
catalytic cap Ocap- Our aim is to obtain solutions for the velocities and rotation for different inclination angles of 
a particle with respect to the wall as a function of 6 and Ocap, and from these solutions identify trajectories that 
can correspond to swimmers being repelled form the surface, skimming over surfaces and stationary in the vicinity 
of the surface. These regimes of motion are especially important in applications, where swimmers are usually in 
the vicinity of boundaries, and engineered motions such as skimming can be useful to precisely direct the swimming 
motion. Tasinkevych et al [77] first examined diffusiophoretic-driven motion of a catalytically-active Janus colloid 
near a wall using the Anderson “slip velocity” framework and constant production of solute for arbitrary active areas, 
wall separation distances and inclination angles. They used boundary element methods and the Lorentz reciprocal 
theorem to obtain particle trajectories, and find that diffusiophoretic-driven Janus colloids can skim along the surface 
for large active areas, and become stationary at very large areas. Our aim in this study is to develope a more complete 
detailed analytical solution using bispherical coordinates, and obtain master curves for the translational and angular 
velocities as a function of separation distance, inclination angle and size of the active area so that trajectories from any 
initial configuration can be computed. Furthermore, we construct a phase diagram to detail the different trajectories 
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FIG. 3. Sketch of a bispherical coordinates 


as a function of the initial orientation and size of the active surface area. This paper is organized as the following, in 
Sec. II we formulate the mass balance conservation for solute and the hydrodynamics problems with a detail analytical 
solution and in Sec. Ill we provide the discussion of the concentration field, swimming velocity and trajectories of 
catalytic colloidal particle close to the wall. 


II. FORMULATION AND ANALYTICAL SOLUTION 


To formulate the problem we begin with the non-dimensional Stokes and continuity equations, 

( 3 ) 

( 4 ) 


V-v = 0, 

V^v - Vp = 0, 


where velocity and pressure are non-dimensionalized with Voo and 


^oo 

R 


respectively, where R is particle radius and 


Voo is defined in Eq. 2 . The sphere and planar wall geometries can be easily described by utilizing a set of orthogonal 
curvilinear bispherical coordinates (a, /3, 0). Bispherical coordinates are related to cylindrical coordinates through the 
following transformations 


e sin a 

cosh P — cos a ’ 
e sinh/3 
cosh P — cos a ’ 


( 5 ) 

(6) 


The range of coordinates are limited toO^a^Tr, —oc</d<oo,O^0^27r, and e is a positive scale factor. 

P = Po represents a sphere whose center is located on the 2 ;-axis at 2 ; = e coth/3o, with radius . ^ ^ . /3 = 0 is a 

smh Po 

sphere of infinite radius and also the plane 2 ; = 0, which is a wall in the present problem. By combining Eq. 5 and 
6 and eliminating /3, it follows that a = constant describes circular arcs which pass through the points z = ±e and 

TT 

are symmetric about 2 ;-axis. The rotation of the surface of a = — around the z-axis produces a sphere with radius e, 

and for the cases of a > ^ and < f spindle-like and apple-like shapes are obtained, respectively. The value a = tt 
corresponds to the line on the z-axis from —e to +e, and the value a = 0 corresponds to the rest of the z-axis (Eig.3). 
Since 2 : < 0 is irrelevant in the physical problem, P just takes non-negative values. 
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The non-dimensional boundary conditions are; 

V|;3=0 = 0, 

=Vs+U-Ff2x(r- rp), 
V —> 0, as —> oo. 


( 7 ) 

(8) 
( 9 ) 


where r is position vector on the particle surface and Vp is the location of particle center of mass which one can show 

6 R 

^0 = cosh“^(—-—) . The main objective is to find translational and angular velocity vector (U, 11) from the fact 
R 

that the particle is freely suspended in a viscous medium. In order to obtain these unknown velocities we need to find 
the distribution of solute around the particle, since the slip velocity at the particle surface is where C is 

NqR 

the solute distribution at the particle surface, non-dimensionalized with ^ . 


A. Concentration Field 

The mass transfer relation which governs the solute distribution in steady state around the colloid is the Laplace 
equation 


= 0 ( 10 ) 

We can easily satisfy the boundary conditions at the wall as well as the particle surface in bispherical coordinates 
((a,/3, 0), where a general 3-D solution of Laplace’s equation can be expressed as [78]: 

C{a,f3,(f)) = v^cosh)d — cos a 

^ ^ 1 _ 1 ( 11 ) 

X y] y] [An,m Sinh(n -)/3 Bn,m cosh(n -)/3] P™(cos a) cos(m0 7 ^). 

m=0 n=m 


Here, P^{cosa) is an associated Legendre polynomial of the first kind, and the constants An^m^ Bn^m and ^rn can be 
obtained upon imposing appropriate boundary conditions at solid surfaces. Using constant flux production of solute 
at the active side of the particle and zero flux at the passive part of the particle as well as the wall surface, the 
boundary conditions can be written in nondimensional form as: 


dc 



( 12 ) 


(cosh pQ — cos a) dC 
€ djS 


/3=/3o 




and /(a, 0) is the coverage function defined as: 



r G Active side 
r G Passive side 


(13) 


(14) 


The relevant recursive relations to find concentration field can be found in the Appendix A, where the conditions 
(12-13) are used to find the unknown coefficients. 


B. Stokes Flow Solution for the Velocity Field 

The number of independent components of the force and torque on the particle is reduced by symmetry to just 
three. The particle has a symmetry axis, h, running through the center of the cap towards the center of the particle, 
which by an appropriate choice of coordinate system can be chosen to lie in the x-z plane, where the z axis is normal 
to the solid boundary plane. In this coordinate system, the particle and concentration fields are symmetric about 





(a) 
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(b) 




FIG. 4. Decomposition of the velocity field into (a) surface diffusiophoretic slip with no net movement (translation or rotation) 
of the particle, (b) translation in x direction with no slip or rotation or motion in the z direction, (c) translation in z direction 
with no slip and rotation and motion in the x direction and (d) rotation with no slip or translation. 


the x-z plane, and there is no force in the ^-direction. Likewise, in the absence of any tendency to rotate out of this 
plane, the angular velocity and torque must be in the ^-direction. Thus 

Fy=T, = T, = 0. (15) 

An alternative argument is to note that the only vectors in the problem are n and z, so the force must be a linear 
combination of these, and the only pseudo vector available is n x z, so the torque must lie in that direction. Choosing 
n as above, the same non-zero force and torque components are obtained. 

The fluid velocity field outside the thin interaction layer is governed by the Stokes flow and continuity equations. 
Due to the linearity of Stokes equation and the boundary conditions associate with it, the hydrodynamic problem can 
be decomposed into four separate problems as illustrated in Fig. 4 

1. “a”: The colloidal particle is stationary, i.e. (14 = 0 and 14 = 0 and ^1^=0), but the fluid slips with the 
diffusiophoretic slip velocity v^. 

2. “b”: The colloidal particle is translating with velocity 44 in the parallel direction along the wall with no rotation 
{fly = 0) or translation in the 2 ) direction (14 = 0)- In this case, there are force and torque shown as Fj and 
Ty exerted on the particle. 

3. “c”: The colloidal particle is translating with velocity 14 in the normal direction to the wall with no rotation 
{Cly = 0) or translational motion along the wall (14 = 0)- In this case, the motion is axisymmetric and there is 
a force on the particle in the ^ direction. 

4. “d”: The colloidal particle is rotating around the y axis with angular velocity fly and there is no translational 
motion in any direction (14 = 0 nnd T4 = 0). In this case, there are force and torque as and T^ on the 
particle. 


The general procedure for solving each one of these problems is the same except for the boundary conditions on the 
colloid surface. In addition, the solutions to problems (h) - (d) have already been examined in the literature, and we 
use the previous results on these problems to validate our calculations [78-82]. 

It was shown previously the slow motion of a sphere in a vicinity of a planar wall can be solved analytically in 
cylindrical coordinate where the components of the velocity can be written as [78]: 




\/cosh/3 — cosa 
e 


T T (^n,m sinh(n+-)/? +TO cosh(n+-)/?) Py (cos a) cos(to</> + 7to), 

m=0 n=m 


X 


(16) 
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sin a 

Vr = - 

2 V cosh 13 — cos a 

oo oo ^ ^ 

X ^ ^ {An,m sinh(n + -)/3 + Bn,m cosh(n + -)/3) P™(cosq;) cos(to0 + 7 ^) 

m=0 n=m 

00 ^ ^ 

+ v^cosh^ - cos a {En,o sinh(n + -)/3 + F„,o cosh(n + -)/3) P^{cosa) cos( 7 o) 


+ E (£;„,„ sinh(n + -)/3 + cosh(n + -)/3) P™+^(cosq;) cos(m</) + 7 ^) 

m=l n=m+l 

^ 00 00 ^ ^ 

+ ^E E (Gn,m sinh(n + — )/5 + cosh(n + — )^) ^(cosof) cos(7ti0 + ' 7 ^)], 

m=l n=:m—1 


(17) 


= \J cosh ^ - cos a (G„,o sinh(n + -)/3 + Hnfi cosh(n + -)/3) P^cosa) sin( 7 o) 

z z 

n=l 

^ 00 00 ^ ^ 

+ ^E E (P„,mSinh(n+-)/3 + P„,mCosh(n+ -)^) P™+^(cosa) sin(m</) + 7 m) ( 18 ) 

m=l n=m+l 

^ oo (X) ^ ^ 

4E E (G„,mSinh(n+-)^ + P„,TOCOsh(n+-)^) P™ ^(cosa) sin(m0 + 7m)], 

m=l n=m—1 


sinh (3 

^ 2 v^cosh P — cos a 

oo oo ^ ^ 

X ^ ^ (^n,m sinh(n+ -)^ + P„,m cosh(n + -)^) P™(cosa) cos(m0 + 7 m) 

m=0 n=m (19) 

+ \/cosh [3 — COS a 

oo oo ^ ^ 

X ^ ^ (C„,mSinh(n+-)/3 + P„,mCosh(n+-)/3) P™(cosQ!) cos(m(/> + 7 ^), 

m=0 n=m 

where An^m^ •••5 Hn^rn ^re 8 sets of unknown coefficients that have to be determined by applying the relevant 

boundary conditions and the equation of continuity. The detailed analysis for problems (a) - (d) are all given in the 
Appendices B-E. 


C. Calculation of hydrodynamic drag forces and torques 

To determine the swimming translational and angular velocities of the colloid, we first take advantage of the the 
fact that the net force and torque exerted on the particle is zero. The drag forces and torques associate with each one 
of the problems (a) - (d) of Fig. 4 are summed and set equal to zero. Having determined the unknown coefficients 
of Stokes equation, we proceed by calculating the force and torque exerted on the colloid due to the diffusiophoretic 
motion of fluid around stationary colloid where, by definition we have, 

F = ^ {a. e„) dS, (20) 

Sp 

where Gn is a unit vector normal to the surface of colloid Sp, and a is the stress tensor 

a = -pi + fi [Vv + (Vv)^]. (21) 

In the axisymmetric problem where the tilt angle E = 0° or 180°, it is obvious by symmetry that the fluid exerts a 
force on the colloid solely in ^-direction and the vertical force exerted on the colloid in terms of coefficients of Stokes 
solution governed by [78] 

CX) ^ 

= ‘I'spl'K e /i ^ ^ (^n,oT-^n,o) {ti —) (A^^o T 

n=0 


( 22 ) 
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where superscript P stands for propulsion. We note here for a Janus colloid in an arbitrary configuration, the vertical 
force on the particle in the ^-direction is only a contribution from m = 0 in the previous formula (16-19). The 
colloidal particle also experiences a force in x-direction which is a contribution from m = 1 component of coefficients 
of equations (16-19). In this case, it can be shown that [58], 

CX) 

= —V^TT e jl Voo (Gn,l+iJn,l) ~ ^(^ + 1) (^n,l + ^np)? (23) 

n=0 

and lastly we need to evaluate the torque that has been exerted on the colloid due to the fluid slip. The general 
expression for torque is: 


T = ® (r - Fp) X (cr. en) dS. 


(24) 


Here, the only non-zero component of the torque is in ^-direction, and it is solely a contribution of m = 1 terms in 
coefficients of pressure and velocity field equations (16-19) and can be written as [58]: 

Ty = V2 71 fi Voo 

oo ^25) 

—n{n + l)(2((7n,i + T^n,i) + cothydo(H^^i -j- H^^i)) — (2n -hi — cothydo)(Gn,i + -f^n,i)- 

n=0 

To calculate the swimming velocity of the particle we use the fact that the colloid is force and torque free. In this 
case, the swimming velocity in z direction can be obtained from 

Ff + pP = 0, (26) 


where the superscripts P refer to normal motion and propulsion problems respectively as noted before. Similarly, 
the translational velocity in the x-direction, Vx^ and the angular velocity in the y direction, 11^, satisfy 

fFj+F^ + Ff = 0, 

\t^ + t^ + tp = o, 

where T, R refer to translation and rotation problems respectively. Expressions for the various forces and torques 
related to problems {b)-{d) can be found in the Appendices C-E. 


D. Reynolds Reciprocal Theorem 

In this section we utilize an alternative approach to determine the diffusiophoretic force and torque on the colloidal 
particle associated with the first problem shown in Eig. 4 . This approach is most useful for cases where the details 
of the flow field are not important, and is based on the Reynolds Reciprocal Theorem (RRT) [83]. Suppose (v',cr') 
and (v",cr") are the solutions of Stokes flow for incompressible fluid within an arbitrary fluid volume E. The RRT 
states that: 


• cr' • v"dS = • cr" • v'dS. (28) 

as as 

Here, the boundary 5E consists of the particle surface, the wall z = 0, and a hemisphere at large distances from the 
particle. Since for a colloidal swimmer the flow field decays as a dipole v ^ and cr ^ r~^ at large distances, the 
contribution from the outer surface integral is always zero and because of the no-slip condition the wall contribution 
vanishes, and hence we are left with an integral on the colloid surface. 

Eor a general configuration where the colloid moves in an arbitrary direction with respect to the wall, similar 
to previous section we need to obtain the propulsion forces, and in z (normal) and x (parallel) directions 
respectively as well as the torque T^, around the y axis. To evaluate these unknowns, we take advantage the 
detailed solution of the three classical problems explained extensively in literature [78-82] and are recapitulated in 
the Appendices C-E. The propulsion force in the z direction can be found by letting (v", cr") be the solution of 
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translational motion of a colloid normal to the planar wall governed by Brenner and Maude (discussed in Appendix 
D) as [79, 80], 


p ^ coshjS - cosadC 

dl) ^/3 ’ ^ct/3 ^ ■ 


da 


Sp 

1, r esinof dC 

/o Jo (cosh pQ — cos a) da 


da. 


/3=/3o 


(29) 


where and are the unit normals in bispherical coordinate and is the stress tensor associate with the problem 
discussed in Appendix D. Inasmuch as is related to an axisymmetric problem, it is independent of azimuthal 
angle and the only term inside the integrand which is a function of 0 is Utilizing the orthogonality, we realize 
that the only surviving term is for the case where m = 0 and therefore the double sum in 11 reduces to one sum. 


277 


da 


d(j) 

^=/3o 


sinof 

2 ycosh /3o — cos a 


oo 

cosh(n + 0.5)yd(A^o cosm^ + Bno sinm0)P^(cos a), 

n=0 


oo 

— sinay^coslpdo^^^o^ cosh(n + 0.5)yd(A^o cosm0 + sinm0) 

n=0 

at this point, we can no longer proceed with analytical calculation and the integral in Eq. 29 needed to be evaluated 
numerically. To do so, we first calculate non-dimensional concentration field and its derivatives as well as surface stress 
tensor component according to appropriate relations and thereafter we undertake integration utilizing trapezoidal rule 
to have the propulsion force in z direction. 

To determine the propulsion force parallel to the wall (in x direction), we let (v", a") in Eq. 28 be the solution of 
translational motion of a non-rotating colloid parallel to a planar wall which satisfy the no slip boundary condition 
on particle surface and the wall. This problem was governed first by O’Neill [82]. Having used RRT, the propulsion 
force in the x direction can be written as. 


dP^ {cos a) 
d cos a 


(30) 


P = - ® • <72 • VCdS = -(/i + /2 + h), 


(31) 


where, 


/l — TT / ((J2 • Ofi • e^-) 

Jo 


€ sin Q( (cosh Pq cos a — 1) 
(cosh Po — cos a)‘^ 


sinof 


COsh(n -h 0.5)/doPnl^n(cOSQ() 


Zy^coshydo — COSQ( 

oo dP^ f cos od 

sma\/ cosh Pq — cos a > cosh(n + 0.5)ydoPni—r- -]da^ 


(32) 


n=l 


h = TT 


[ ^ CQsh Po — cos a cosh.{n P O.b)P qBP^ {cos a)]da^ 

Jo (cosh Po — cos Q() 


(33) 


^ . e sinh/dosin^Q( . sin a ^ 

h = -7T ((72 • On • e^)--—- ;2[ o / n /Q 

Jo (cosh Po — cos a) 2v cosh po — cos a 

oo dP^ f cos od 

— sin Q(V cosh/do — cos a > cosh(n + 0.5)/doPni —t - -]da, 

^ d cos a 


cosh(n + 0.5)ydoPni^n(cos a) 


n=l 


(34) 


where <72 accounts for the stress tensor associate with the O’Neill’s analysis [82] (see Appendix C) and B^i can be 
found according to the recursive relations given in Eq. 37 and 38. 
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Lastly, to obtain the propulsion torque around the y axis, we let (v", a") in Eq. 28 be the solution of rotation of 
a non-translating colloid around the y axis in a proximity of a wall. This problem was derived in detail by Dean and 
O’Neill [81] (see Appendix E). Having used RRT, the diffusiophoretic torque around the y axis can be written likewise 
as. 


T^^ = -e(/i+J2 + /3), (35) 

where Ji, I 2 and Is are given in 32 - 34 except the stress tensor (73 associate with the Dean and O’Neill’s solution 
should be replaced [81]. 

Having determined the propulsion force and torque, we again utilize the fact that the total force and torque on colloid 
should be zero and via a force/torque balance in three directions, one can obtain the swimming velocity and angular 
velocity of the colloid by using Eq. 26 and 27. Our results from RRT analysis are in complete agreement (variations are 
less than 1 % in all separation distances and coverages) with the direct solution of Stokes flow given in previous section. 


III. RESULT AND DISCUSSION 

A. Concentration distribution 


We begin with the case of an inclination angle of 0° or 180°, where the concentration field around the colloid is 
independent of azimuthal angle (j). When the active section of the colloid faces away from the wall, E = 180°, the solute 
distribution does not change drastically as the colloid approaches the wall since the solute can diffuse freely away from 
the colloid to distant regions far from it. On the contrary, if the active side of the colloid faces the wall, 5 = 0°, as 
the particle approaches the wall the concentration field is highly distorted by the presence of an impenetrable barrier 
resulting in a very high solute concentration in a thin gap region adjacent to the wall (see Eig. 5 for the case Ocap = 90 
)• 



Eor all other colloid orientations, 0° < 5 < 180°, the concentration field is fully three dimensional in bispherical 
coordinates C{a^P^(j)). Representative concentration fields in x-z plane around a Janus particle {Ocap = 90), for 
various inclination angles and a fixed separation distance J = 1, are shown in Eig. 6 . The effect of the waii is 
again more pronounced when the active sections faces the waii and as the inciination angie approaches zero, where 
the soiute concentration is augmented in the thin gap. In generai, the concentration fieids are no ionger symmetric 
around the coiioid axis (see Eig. 2) which resuits in a rotation about it which we wiii discuss in detaii presentiy. 
The asymmetric distribution around the coiioid axis is more pronounced as the inciination angie approaches to 5 = 90°. 

The degree of concentration asymmetry about the coiioid axis can be aitered as weii by varying the extent of the 
active section of the coiioid surface. Eig. 7 iiiustrates the concentration fieids in the x — z piane for the case of smaii 
active coverage, Ocap = 40°, as weii as for a greater active section, Ocap = 140°. Note that for Ocap = 140°, there is 
a substantiai increase in soiute concentration in the gap reiative to the case where Ocap = 90° for simiiar inciination 
angies. Thus, we see that the effect of soiid waii for coiioids with high active coverage is more significant compared 
to the ones with smaiier coverage. 
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FIG. 6. Non-dimensional concentration field for a Janus particle (Ocap — 90) in the x-z plane for different inclination angles 
(S) when the particle is one radius away from the wall. 


B. Swimming velocities and trajectories 


Swimming normal to the solid wall 


Fig. 8(a) illustrates normalized swimming velocities (non-dimensionalized with the swimming velocity Voo for an 
infinite medium) as a function of non-dimensional separation distance 6 (non-dimensionalized with the colloid particle 
radius i?), for different inclination angles for a half catalytically active Janus particle {Ocap = 90°) . When the initial 
orientation is axisymmetric, S = 0° or S = 180°, the particle translates in the ^-direction without rotation, but the 
behavior differs in the two cases. When the active side of the colloid faces the wall the solute concentration builds 
up in the gap region, resulting in high gradient and a large slip velocity, leading to a large diffusiophoretic propulsive 

w 

force driving the particle away at an enhanced velocity > 1. Of course at large separations the wall lose relevance 

Foo 

and the velocity relaxes to an infinite medium value = 1. In the opposite case when the active side faces away 

Fcx3 

from the wall, the concentration filed and diffusiophoretic force is little altered by the boundary but the hydrodynamic 
lubrication force [80], growing as ^ dominates at small gaps and brings the particle nearly to rest. Again, at large 

w 

distance from the wall, the infinite medium result is recovered, = — 1 in this case. These general remarks still 

Foo 

apply for other values of Ocap except for the case of very high coverage which will be discussed later. 


2. Repulsion and refleetion of a eolloid from wall 


When the orientation of particle is no longer axisymmetric, the concentration field is fully three dimensional, so 
that the particle moves in the x-z plane and rotates about the ^-axis. The wall is effective only at separations 6 < b 

Vx Vz 

and at larger distances the infinite medium result applies with (t^, tw) ^ (sinS, cosS). 

FcxD Fqo 

Consider first the case of a hemispherical active area {Ocap = 90°). The trajectory of a Janus colloid initially 
located at J = 1.0 with an inclination angle of S = 70° is plotted in figure 8(d) where the colloid is repelled from 
the wall without considerable rotation. This trajectory is a typical of all partially active catalytically colloids with 
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id) 2 = 0 = 140“ {e) 2 = 90 ,0^„^ = 14O" (/) 2 = 180“,0^.^^ = 140“ 

FIG. 7. Non-dimensional concentration field for low active coverage {Ocap = 40) and high active coverage (Ocap = 140) in the 
x-z plane for different tilt angles (S) at fixed distance from wall (5=1. 


different surface coverage and initial separations comparable to the particle radius with initial orientation S < 90°. 
This behavior can be easily understood from the results given in Fig. 8, where we see that for E < 90°, Vz > 0 for all 
separation distances and therefore the colloid would never approach the wall. 

For initial orientation E > 90°, the diffusiophoresis force is towards the wall and the behavior is more interesting. 
Figure 8(e) exhibits two trajectories of a Janus colloids Ocap = 90° with a same initial separation from the wall, 
(5 = 3, but different inclination angles S = 120° and S = 150°. In both cases the colloid approaches the wall with a 
negligible counterclockwise rotation, followed by a fast counterclockwise rotation in the vicinity of wall which causes 
the active side of the colloid to face the wall causing it to be eventually repelled from the wall. As the initial inclination 
angle approaches 180°, the colloid experiences reflection from the wall at smaller separations. Note that the angular 
velocity grows noticeably as the colloid approaches the wall (see Fig. 8(c)). The final configuration of the particle 
has inclination angles less than but close to 90° which causes the colloid to move away very slowly in the positive ^ 
direction. Finally we ask whether there is any other trajectory beside reflection from the wall for a Janus particle 
{^cap = 90°). We first note that at each separation distance from the wall, there is a specific inclination angle for which 
Vz = 0. Furthermore, since ^4 < 0 at all separation distances, the colloid rotates counterclockwise, which decreases S, 
which leaves the colloid at an orientation where Vz > 0 and it tends to move away. 

To explore the other possible scenarios for the swimmer against planar wall, we vary the active surface coverage of 
the colloid. Fig. 9 illustrates non-dimensional swimming velocities as a function of normalized separation distance for 
various inclination angles for a colloid which is mostly passive; Ocap = 40°. From these curves, it can be concluded 
again that repulsive interaction of solute and particle eventually brings about a net repulsion of the colloid from the 
wall. A typical trajectory of a particle with low active coverage Ocap = 40°, starting from separation distance S = 3.0 
and inclination angle E = 140°, is shown in Fig. 9(d). We see that the colloid approaches the wall with a small 
counterclockwise rotation, but due to the strong hydrodynamic repulsion in the vicinity of the wall, the colloid rotates 
into a configuration where Vz becomes positive and the particle ultimately is reflected form the wall. According to 
the angular velocity curves in Fig. 9(c), an orientation with positive rotation is possible but since at that orientation 
the velocity in the 2 ; direction is positive the colloid moves away from the wall and subsequently its rotation decays 
to zero. Our analysis suggests that repulsion from the wall is typical for colloids with initial orientations of S < 90° 
for any active coverage section Ocap: reflection from the wall is obtained for orientation angle E > 90° for active 

coverage of 0° < Ocap < 115°. 
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FIG. 8. 8(a) Non-dimensional swimming velocity in the z direction, 8(b) non-dimensional swimming velocity in the x direction, 
8(b) Non-dimensional angular velocity as a function of separation distance S for half covered Janus particle (Ocap — 90°). 
Sample trajectories for half covered Janus particle {Ocap = 90°) 8(d) initial conhguration: J = 1.0, S = 70° and 8(e) initial 
configuration: S = 3.0, S = 120°, S = 150°. 


3. Skimming of a colloid along planar wall 


As the coverage increases above Ocap = 115° other types of trajectory appear for the colloids with the active side 
facing away from the wall. In Fig. 10(a)-10(c) we show the velocities for the case Ocap = 120°, indicating that there 
is always a specific configuration of colloid (separation distance and inclination angle, 6s and S^) in which Vz and 
fly are close to zero, so that the colloid slides along the wall while maintaining its specific orientation and separation 
from the wall. If the particle initially has an inclination angle greater than 90°, eventually it slides along the wall 
with trajectories similar to the ones shown in Fig. 10(d). This figure indicates different initial configurations leading 
to identical sliding behavior. Generally, having known the initial configuration of colloid, the final orientation of 
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FIG. 9. 9(a) Non-dimensional swimming velocity in the z direction, 9(b) non-dimensional swimming velocity in the x direction, 
9(b) Non-dimensional angular velocity as a function of separation distance 5 for low coverage colloid (Ocap — 40°) and 9(d) 
sample trajectories for colloid with coverage {Ocap = 40°) and initial conhguration: 6 = 3.0, S = 140°. 


the colloid with respect to the wall can be extracted from the curves in Fig. 10(e)-10(f), which give the variation of 
skimming distance and tilt angle with surface coverage. In this steady condition, the inclination angle and separation 
distance remains constant while propulsion force in the 2 : direction and propulsion torque are close to zero. As the 
area of active cap increases, particle skims in the smaller separation distance and inclination angle increases and 
approaches 180°. 


4- Stationary Particle 

We have seen that when the surface coverage is not too high, Ocap < 145° and the inclination angle is close to 
180°, i.e. the active side facing away from the wall, the solute repulsion drives the particle towards the wall at all 
separation distances. Surprisingly, when the coverage exceeds this value of Ocap the particle is instead repelled at 
small separations! The reason is that if the particle is near the wall the active surface region extends into the thin 
gap region and produces a local repulsion which reduces the total diffusiophoretic force to values comparable to the 
lubrication resistance. The concentration field in such case is given in Fig. 11(a). When a particle of high coverage 
approaches a wall at a large value of S, this effect produces a positive rotation (Fig. 11 (d)) which drives the orientation 
to S = 180°, where the force in the x-direction vanishes and the particle moves to rest (Fig. ll(b)-Fig. 11(d)). Two 
examples of such stopping trajectories are given in Fig. 11(e). We see in these figures at I 4 = 0 since E = 180°, but 
Vz = 0 only at a certain value of S, so the final separation distance Sf depends on Ocap^ as indicated in Fig. 11(f), but 
as the trajectory plot indicates Sf is insensitive to initial orientation. 
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(c) 





FIG. 10. 10(a) Non-dimensional swimming velocity in the ^-direction, 10(b) non-dimensional swimming velocity in the X- 
direction, 10(c) Non-dimensional angular velocity as a function of separation distance 6 for high coverage colloid {Ocap = 120°), 
10(d) sample trajectories for colloid with coverage (Ocap — 120°) and initial configuration: (5 = 1.0, S = 120° and 5 — 4.0, 
S = 150°, 10(e)-10(f) equilibrium skimming distance and tilt angles respectively as a function of surface coverage. 


IV. SUMMARY AND CONCLUSION 

This study has aimed to study the self-diffusiophoretic motion of a catalytically active swimmers near a solid planar 
boundary. In self-diffusiophoresis, the colloid creates and sustains a solute gradient as a means of autonomous motion, 
and the colloid functions as an engine or motor. The hydrodynamic and mass transfer equations which describe the 
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FIG. 11. 11(a), concentration field for a final configuration of particle with {Ocap — 150°), 11(b) Non-dimensional swimming 
velocity in the ^-direction, 11(c) non-dimensional swimming velocity in the X-direction, 11(d) Non-dimensional angular velocity 
as a function of separation distance S for very high surface coverage of particle {Ocap = 150°), 10(d) sample trajectories for 
colloid with coverage {Ocap = 150°) and initial configurations: (5 = 1.0 , S = 120° and S = 4.0 , S = 140° and 11(f) equilibrium 
fix position distance of particle to the wall as a function of coverage. 
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self-phoretic motion of the colloid are formulated utilizing slip velocity approach [46] in the small Damkohler Da 1 
and Peclet Pe 1 number regimes. The hydrodynamic problem was solved exactly, both by directly solution of the 
boundary value problem and via Reynolds Reciprocal Theorem. 

We have shown the complex dynamics of catalytically active swimmers near a solid boundary is controlled by the 
orientation of the colloid relative to the solid surface, E as well as the active coverage of the catalytic cap, Ocap- For 
situations where the tilt angle equals to 0° or 180°, the motion is axisymmetric and simple rectilinear motion in the z 
direction results, where symmetry insures that the colloid does not rotate and its trajectories are perpendicular to the 
wall. When the active area faces the wall, the solute concentration rises in the gap between the wall and the colloid 
due to the zero flux condition at the wall surface, and the particle experiences a boost in diffusiophoresis relative to 
the propulsion in an infinite medium. This effect is opposed by viscous drag, as modified by the presence of the wall, 
and the balance between these two effects determines the swimming velocity. 

When the active area of the swimmer is not aligned perfectly towards or away from the wall, the inclination generates 
a diffusiophoretic propulsion along the wall, perpendicular to the wall as well as a torque perpendicular to the plane 
of symmetry passing through the particle. For 0° < S < 90° and for all coverages, the net diffusiophoretic propulsive 
force is in the positive z direction, since the orientation is such that solute accumulates between the wall and the 
particle, pushing the particle away from the wall with small rotation that vanishes as the particle moves further from 
the wall. For E > 90°, the diffusiophoretic propulsion in the z direction may vanish at a certain critical tilt. In this 
case, the reactive cap is obliquely oriented to the surface so that the accumulation in the gap is reduced and the 
increased concentration above the colloid balances the reduced upward propulsion from the bottom hemisphere. In 
this state the particle, at least instantly, skims along the surface. For a given Ocap, the critical tilt will be a function 
of the separation distance, and this relationship will vary with the cap size. We can anticipate that as the cap size 
increases, the critical angle for a fixed S will increase as more of the cap is required to face away from the wall to 
reduce the diffusiophoretic boost due to the accumulation in the gap between the wall and the particle. For particles 
with low to medium surface coverages at that critical tilt angles which is a function of separation distance at all values 
of 6 rotation is counterclockwise at the critical tilt angles, which leads to repulsion of the particle from the wall, along 
with a reduction of the inclination angle. For particles with high coverage there is a unique configuration at each Ocap 
at which the normal force and torque are close to zero, so that the particle skims at that steady configuration. Finally 
particles with very high Ocap attain their fixed position with an axisymmetric orientation and the particle will halt. 

Based on the results we have obtained, we can construct a phase diagram which gives the long-time behavior of 
a partially active catalytic particle near a solid planar wall: Fig. 12. We believe the phase diagram can have an 
important implications on designing the next generation of micro-nano engines for various biological applications. For 
example, the skimming regime is an area of great interest since it offers the possibility of steering a particle along 
a solid wall from one place to another. Likewise, the stationary regime has obvious applications to drug delivery at 
particular positions. 

One important issue not considered here is the effects of thermal (Brownian) fluctuations on colloidal particle 
motion. In general, one would anticipate that large colloidal particles would be resistive to thermal fluctuations 
and the current results would apply, while smaller particles might behave differently, in particular the existence of 
skimming or stationary states may be problematic for sufficiently smaller colloids. We intend to return to this issue 
in a future publication. 



20 



FIG. 12. Types of particle trajectory as a function of tilt angle and active cap size. 


V. APPENDIX 

We first discuss the concentration field calculations, and then the strategy used to obtain solutions to the problems 
of Fig. 4. As discussed before, due to the linearity of Stokes equation and boundary conditions, it is possible 
to decompose the general problem into some subproblems and by superposing the solutions to the problem we can 
construct the flow associate with the general problem. The general solution for the flow field around a spherical 
particle close to a wall was proposed by Leal [78], given in Eq. 16 - 19. In the following, we review the boundary 
conditions associated with each problem shown in Fig. 4 and provide the analysis and related recursive relations. 


A. Concentration field around the particle 

The boundary condition (12) is satisfied trivially when 

— b ^5 ^ ^ b* 


(36) 


Furthermore, we choose the axis in such a way that the active cap is completely symmetric about the x-z plane 
and therefore 7^ = 0. By imposing boundary conditions (13 -14) at the colloid surface, and utilizing the appropriate 
orthogonality and recursion relations among the associated Legendre polynomials, we have 


1 — 1 3 — 

[-nsinh(n - -)/3o Bn-i,o + (n sinh(n - -)(3o + (n + 1 ) sinh(n + -)/3o)Bn,o 
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(2n+l) 
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sin a (cos a) da d(j) (m. = 0), 


(37) 


1 — 1 3 — 

[-(n - m) sinh(n - -)/3o + (n sinh(n - -)/3o + (n + 1) sinh(n + -)(io)Bn,r. 
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( 38 ) 
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The special cases S = 0^ and E = 180^ are two axisymmetric situations (where the active side faces toward and away 
from the wall, respectively) where the concentration is independent of azimuthal angle and the only relevant term 
is m = 0. 

The infinite series solution is approximated by truncating the series for large values of n and m, based on the 
fact that the coefficients vanishes for sufficiently large number of terms, and the remaining unknown coefficients are 
obtained by solving equations 37 - 38. Lastly, we note that more terms in the truncated series are required for small 
separation distances to reach numerically accurate and consistent result for the concentration field. 


B. Stationary sphere near a wall with slip on the sphere surface 


A spherical particle is located in a vicinity of an infinite plane wall z = 0, where the sphere is stationary but the 
fluid slips with velocity on its surface. In this case non-dimensional boundary conditions are, 

= 0, (39) 

= -VsC, (40) 

V —y 0 , (25 (41) 

The no-slip boundary condition at the wall can be expressed in cylindrical coordinate as. 


1/3=0 ~ '^01/3=0 ~ '^2; 1/3=0 — 0- (42) 

We proceed first by determining the velocity field for m = 0, which corresponds to the axisymmetric problem, 
by evaluating An,o, 5n,o, •••, ^n,o- The slip velocity at the particle surface is proportional to surface gradient of 
concentration VgC, where Vg = I — nn. Using the relations between tangent unit vectors to the surface of sphere in 
bispherical coordinate and their counterparts in cylindrical coordinates one can show that 


l/3=/3o 


^ ^ r (cosh 3o cos a — 1) sin a ^ ~ ^ ^ 

= - VgC-er = - p - y] ^ri,o cosh(n+ -)^o Pn{cosa) 

2e A/cosh Po — cos a ^ 2 
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(cosh/3o cos a - 1) i/cosh/3o - cosa ^ ~ Ino ni/ m 

-2^-B„,o cosh(n + -)/3o P„(cosq;)J, 
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(43) 


Wl/3=/3o ~ • 60 — O 5 


(44) 


l/3o 


^ . sinh p sin^ a ^ \ 

sC e^= -— / , ^ ^ cosh(n-h -)do Pn(cos(a) 


2e \/cosli)3^^^cosa 

sinh/3o sin a \/cosh/3o - cosa ^ ~ 1 ^ 

H- - - 2 ^ Bnfi cosh(n + -)/?o -P„(cosq!)J. 

Axisymmetry implies that 70 = 0, and 


n=l 


(45) 


^n,0 — P-n,0 — 0* 

The no slip condition at the wall (42) yields the following recursive relation 

1 ^ (n-l) 


1 

'2(2n- 1) 


Pn-1,0 + 


2(2n + 3) 


B, 


n-|-l,0 


(2n - 1) 


F IF (^ + 2) p _ „ 

-^n— 1,0 T -^n ,0 — C! 


(2n + 3)' 


72 > 1, 


and 


(46) 


(47) 


-Dn.o = 0 ;n>0. 


(48) 
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Applying the slip velocity condition at the particle surface, Eqs. 43 and 45, to the solution of equation of motion for 
m = 0, Eq. 17 and 19, we obtain 


- 2[2n- 1) + Bn-ifi cosh(n - i)/3o] 

13 3 

+ 2(2n + 3) sinh(n+ -)^o +-B„+i,o cosh(n + -)/3o] 

- (2n + Fn-ifi cosh(n - i)/3o] 

+ cosh/3o[£^„,o sinh(n + ]^)(5q + F„,ocosh(n + i)^o] 

^ 2) 3 3 

- (2n + 3) sinh(n+ -)/3o + F„+i,ocosh(n + -)/3o] 

= -4[cosh/3o cosh(n - 5)/3o ( -B„_ 2 ,o 

+ cosh(n- h/3o ( ^ (1 + 2(n - 1)(1 + cosh^/3o))) Bn-ifi 

- cosh/3o cosh(n + i)/3o + 2) -Bn,o 

-cosh(n+ ^)/3o (—2—(l-2(n + 2)(l + cosh^/3o)))-B„+i,o 
2 2n + 3 

5 —(n + 2) ~ 

+ cosh/3o cosh(n + -)/3o —-^)B„+ 2 .o] (n > 1), 

2 2n + 3 


(49) 


sinh^O [^^,o ^ ^ cosh(n + i)/3o] 

77/1 1 

- (2n - 1) sinh(n - -)^o + -D„-i,o cosh(n - -)^o] 

+ cosh/JoiCn.o sinh(n + i)/3o + -C>n,o cosh(n + ]j)Po\ 

3 3 

- (2n + 3) sinh(n+ -)/3o + L>„+i,o cosh(n + -)/3o] 


sinh^Or . / 3 n (n - 1) ~ 

[cosh(n - -)/3o —j-) Bn- 


2e 


2n ■ 


2,0 


- cosh(n - cosh^o Bn-i,o 

1 1 n(2n^ + 3n-l) (n +1) (-2n^ - n + 2) ~ 

- cosh(n + p/3o (-)) 

+ cosh(n+ cosh^o -Bn+i,o 

- cosh(n + + -Sn+2,0] (n > 0). 


Lastly, the continuity equation gives two more recursive relations 

— 1 5 1 

^- 77 / + 2 ^^,0 ^ 2 ^ ^ n + 1,0 ~ ^ ^ n - 1,0 + ( 277 / + l ) Dn ,0 

— (tT/ + — 77/(77/ — 1)E^_1^0 + 277/(77/ + 1)E^^0 “ (^ + 1)(^ + 2)E^+1^0 = O 5 


— 1 5 1 

— 77/ + 2 + 2 ^n+1,0 ~ ^ Cn-1,0 + (277/ + l)Cn,0 

— (tT/ + l)(7n+l,0 ~ '^(^ ~ + 277/(77/ + 1)E’^^0 “ (^ + 1)(^ + 2)E’^+1^0 = O5 


(50) 


(51) 


( 52 ) 
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For m > 1, no-slip boundary condition at the wall gives two sets of recursion relations 

1 


2(2n- 1) 

{ti -|- TTl -|- 2) 


-^n—l,m “1“ 


1 ^ (n - TO - 1) ^ 

. -1 \ -^n—1,' 


(2n + 3) 

(n — m){n — to + 1) 


2(2n + 3)“”’^"’"‘ (2n-l) 

-Ci+i,m = 0 ;TO>l,n>TO+l, 


2(2n- 1) 
(n — TO + 1) 


(n + to) (n + TO + 1) 
2(2n + 3) 


(fl “h TTL} 

Hn—l,m -^n,m 77^ j y7-^n-|-l,m = 0 f- 


(53) 


(54) 


{2n — 1) (2n + 3) 

The no-penetration condition at the solid surface of wall, Vz\^^q = 0, is trivially satisfied if 

Dn,m = 0 (55) 

Next, we impose the slip boundary condition at the colloid surface (Eq. 40). Transforming the unit vector from bipolar 
spherical to cylindrical coordinates. 


^»'l/3=/3o “ “ 


(cosh/3o COSO! - 1) ^ , 

^ / ('^m T Cm T ^m) COSy^TTlCp -|- 

(cosh po — cos a) 


m=l 


= - VsC • 60 = - y] sin(TO(/) + 7m), 


m=l 

oo 


where 


, ^ ^ sinh Po sin a / n > n / / n 

= - V,C • e, = 2^ + Cm + /^m) COs(toC> + 7m). 

^ ^ ’ 777=1 

sin Of ^cosh Po — cos a 


= 


2e 


nm cosh(n + -)po P^(cos(a), 


Cm — 


(cosh Po — cos a) 2 

^2^ 


Pnm(^ + ^)(^ - ^ + 1) COsh(n -h -)ydo ^COSQf), 


n=m—1 


(coshydn — COSQf) 2 / 

Km. = -- y] Bnm COSh(n + -)/ 3 o (cOS «), 

n=m-|-l 


2e 


(cosh Po — cos a) 2 ~ , 1 ^ ^ ^ 

I'm = -^- 2^ Bnm m cosh(n + -)/3o Pn {cosa). 

—e sinof 2 


(56) 

(57) 

(58) 

(59) 

(60) 

(61) 

(62) 


Inserting those boundary values from Eqs. 56-57 into the solution of the Stokes equations in cylindrical coordinate 
(Eqs. 17-18) and using orthogonality along with the recursion relations between associated Legendre polynomials, two 
sets of recursion relations among unknown coefficients found. Eor m> 1 and I > m + 1, 

- 2(2/- 1) + Pi- 1 ,mCOSh(l - i)^o] 

13 3 

+ „/„r , [^i + l,m sinh(i + -)l3o + P; + i,mCOSh(i + -),5o] 


2{2l + 3) 

(1 — ^ — 1) I I 

( 2 ) 1) siiih(/ 2^^*^ Bi—i^m cosh(^ 2 )^*^) 

+ cosh^o[Pi,m Sinh(l + i)^o + Pi,mcosh(l + i)^o] 

(; + to + 2)^^ ■ ui, , 3,^ , ^ ui, ,3,^1 (2; + i)(;-to-i)! 

“ (2i + 3) + 2*'’“ + + 2>'’“1 = 2(i + m+l)! 

TT 

/ / / cosh Po cos ol 1., . / ~ ~ . 7->777 -i-i / \ • 7 

-(( ' - )[Vm + Cm + ^m) + VCQShpO " COS Of Um) Pi (COS Of) SlUadOf, 

V cosh Po — cos a 


( 63 ) 
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while for m > 1 and I > m — 1, 


(I — m)(l — m1), , , ,, 

-- -[Ai_i^rn sinh(; - -)^o + cosh(Z - -)/3o] 

+ - 2(21 + 3) - ~i^W,rn smh(Z + -)/3o + Bi+i,mCosh(Z + -)/3o] 

- sinh(; - ^)po+Hi-l,mCOSh{l - ^)/3o] 

+ cosh/3o[Gz,m sinh(^ + i)/3o + cosh(^ + i)/3o] 

(^ + to) • v,/; I ^^/9 < ZJ V,/-; I ^^/9 1 _ ™ 

■ + 2^^°^ “ 2(; + to-1)! 


/-« 


cosh Pq cos Of — 1 
v^cosh Pq — cos a 


)('^m + Cm + ^m) ~ V cosh Pq — COS Of z/^) PJ^ ^(cosQf) sin Of da. 


(64) 


Inserting the boundary condition for 14, Eq. 58, into the z component of Stokes equation (Eq. 19) we have. 


[An,m sinh(n + ^)/3o + cosh(n + i)^o] 
in — m) . _ 

- ^2^ _ [Gn-i,m smh(n - -)/3o + cosh(n - -)/3o] 

+ cosh/3o [Cn,Tn sinh(n + i)/3o + £>„,„ cosh(n + i)/3o] 

- ^CCtC-3)^^ [Cn+i,m sinh(n + ^)Po + -D„+i,m cosh(n + ^)/3o\ 

sinh/3op 3 (n - to - 1) (n - to) ~ 

=-5^[c»h(n- jj/So- 

1/ 1 2(n — m) (n — 1) ~ 

- cosh(n- -)/3o cosh/3o -----5n-i,o 

2 2n — 1 

- cosh (..+Imi <"+“>»"+“ ~ “** 

^ (n — m + l){{n — m + 2) — 2n{n + m + 1)) ^ 1 


2ti “h 3 


(2n + l) 


) ^n,0 


, , 3^o uo 2(n + TO+1) (n + 2) 5 

+ cosh(n+ -)/3o cosh/3o -z- Bn+i,o 

2 2n + 3 

5. ^ (n + m + 1) (n + m + 2) ~ 

- cosh(n + -)^o-^-^n+2,0 ;m>l, n > m. 

2 2n + 3 


(65) 


To complete the solution two more relations are needed, and these are obtained from the continuity equation. Using 
Eqs. 17-19 the continuity equation (Eq. 3) can be expressed in the form of a double series of terms of the form 
cosh(n + ^)f3 P^{cosa) cos(m0 + 7^) and sinh(n + ^)f3 P^{cosa) cos(m0 + 7^), setting the coefficients to zero we 
find for m > 1 and n > m 


^ (n TTi) ^ ^ (n + 777/ + 1) (ti tti) + {2ti + 

-{nPmP l)T>n+i,m - ]^{n-m-l){n- m)En-i,m P{nPmPl){n- m)En,m 
~ 2 P mP l){n p m p + -Gn-l,m — Gn,m + -^Gn+l,m = O5 


( 66 ) 
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^ {tI Tfl) H- ^ Bji jYi + ^ {tI -|- Tfl H- 1) Bjij^\ jYi {in ini) Cn—l,m “1“ (271 -|- l)C^n,m 

- (tI + 771 + l)Cn+l,m “ ^ (tI - 711 - l)(71 - m)Fn-l,m + (tI + 711 + l)(71 - m)Fn,m 
~ 2 + 771 + 1)(71 + 771 + 2)Fnni,m + ~ ^n,m + = O 5 


(67) 


We have obtained 8 independent recursion relations for ^-independent and 8 more for ^-dependent parts of the field, 
which are solved simultaneously to find the unknown coefficients An ,0 5 ^n, 0 v 5 ^n,o and ^n,mv 5 ^n,m- 


C. Translation of a non-rotating sphere parallel to the wall 


Problem (h) concerns the translational motion of a non-rotating solid sphere parallel to a solid wall, which has been 
solved by O’Neill [82]. In this case, the boundary conditions on the wall are 


K’1/3=0 ~ ^(t>\j3=0 ~~ "^2; 1/3=0 — O 5 

and the no-slip condition on the particle surface is 


( 68 ) 


Vr\p^= COScf), 

(69) 

y<}>\/3o = -sin0, 

(70) 

7^ 

0 

II 

0 

(71) 


One may easily show that in the general solution of the Stokes equation just the ttt, = 1 term survives and furthermore 
7o = 0. Using Eqs. 69 - 70 we have two recursion relations: 


1 


2(2n-l) [■ 

1 

2(2n + 3) [' 
(n-2) 


An- 1,1 sinh(n - cosh(n - ]^)(5q 

3 3 

An+1,1 sinh(n + -)/3o + Bn+i,i cosh(n + -)/3o 


(2n - 1) + Pn- 1,1 cosh(n - i)/3o 

+ cosh^o En,i sinh(n + i)/3o + Fn,i cosh(n + i)/3o 
- /2n+^3) sinh(n + 5)/3o + -Fn+i.i cosh(n + 


= 0, (n>2), 


(72) 


(n — 1) n 


2(2n - 1) 

(n + l)(n + 2) 


An- 1,1 sinh(n - i)/3o + Bn-i.i cosh(n - i)/3o 


2(2n + 3) 
n 


3 3 

An+1,1 sinh(n + -)/3o + Bn+i,i cosh(n + -)/3o 


(2n - 1) 
+ cosh ^0 
(n + 1) 


Gn- 1,1 sinh(n - i)/3o + Hn-i,i cosh(n - i)/3o 
Gn,i sinh(n + ^)/3o + J7„,i cosh(n + 


2 V 2 


3 3 

( 2 n + 3 ) sinh(n+-)/3o + -f7„+i,icosh(n+-)^o 

^ j ^ + ^ 

277 / — 1 / \ 277 / -|- 3 


cosh/3o e ("+ 2)^0 - 


"(’^+i)/^o 


(tT/ > 0). 


( 73 ) 
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From the boundary condition Eq. 71, the no-slip boundary condition and the no-penetration condition at the wall 
surface, we obtain 


sinh Po 


An,i sinh(n + ^)/3o + Bn,i cosh(n + i)/3o 


_ jn-l) 
(2n - 1) 

+ cosh /3o 

{n + 2) 


Cn- 1,1 sinh(n - i)/3o + cosh(n - i)/3o 

Cn,i sinh(n + i)/3o + D„,i cosh(n + 


(74) 


3 3 

(2n + 3) sinh(n + -)(5q + Z)„+i,i cosh(n + -)/3o 


= 0, (n > 1). 


2(2n- 1) 


B 


(n- 1) n ^ 
2(2n- 1) 


I 1 p (»-2) p 

2(2n + 3)^"+^’^ (2n-l)'^" 

(n + l)(n + 2) ^ u 

2(2n + 3) (2n - 1) "“hi 

+i7„,i - =0 (n > 0), 

(2n + 6) 


(75) 


(76) 


-D„,i =0 (n > 1). (77) 

To complete the solution we must satisfy the equation of continuity, which from Eqs. 66 -67 with m = 1 can be written 
as: 


— 1 5 1 

^-(n — 1) An-i^i + - An^i + -{n -h 2) — (n — 1) Dn-i^i 

-|- (2n + l)Dn^i — (n + 2)D^+iq — -(n — 2)(n — l)E^_iq + (n + 2)(n — (78) 

~ 2 ^ ^ + -Gn-1,1 ~ Gn,i + -Gn+i,i = 0, 

—1 5 1 

^-(n — 1) - 5n,l + 2^^ ^ ^n+1,1 ~ ~ 1) Cn- 1,1 

(2n + l)Cn,i — (n -h 2)(7n+i,i — -(u- — 2)(n — l)F^_iq + (n -h 2)(n — (79) 

~ 2 ^ ^)(^ 3)F^+iq + -Hn-1,1 — Hn,l + -i7n+l,l = 0- 

We obtain the unknown coefficients from simultaneous solution of the above equations. The fluid exerts a force on 
the particle, given by [82]: 


CX) 


— —V^TT e fl Ux {Gn,l-\-Hn,l) — n(n + 1) {An,l + ^77,1)5 

n=0 

where Ux is the particle velocity and similarly the torque exerted on the colloid is in the y direction, given by 
Ty = Ux y] [(2 + (n(n+ 1) (2C„,i + coth/3oA„,i) + (2n + 1 - coth/JoGn.i) 

^ n=0 

(2 - (n(n + 1) {2Dn,i + coth/3o-B„,i) + (2n + 1 - coth/3o-f7n,i) • 


(80) 


( 81 ) 
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D. Translation of a sphere normal to the solid wall 

This problem has been solve by Brenner and Maude separately from stream function solution because of axisymmet- 
ric al nature of problem [79, 80]. Here we address this problem by directly solving the Stokes equation in cylindrical 
coordinates, using the expressions in Eqs. 17 - 19. The no-slip boundary conditions at the wall and on the particle 
are 


K|/3=o ~ "^01/3=0 ~ — 0, 

= r 


(82) 

(83) 

(84) 


and from the above boundary conditions and by symmetry we see that the solution only requires m = 0 and 70 = 0 . 
Using the expression for the generating function of Legendre polynomials, we can expand boundary condition at the 
colloid surface, and after some algebraic manipulation we find 


1 


2 ( 2 n- 1 ) 

1 

2(2n + 3) !_ 
(n- 1 ) 


An-ifi sinh(n - i)/3o + B„_i,ocosh(n - i)/3o 

3 3 

An+ifi sinh(n + -)/3o + cosh(n + -)/3o 


( 2 n - 1 ) L' 

■ cosh /3o 

(n + 2 ) 
(2n + 3) [■ 


-E’n- 1,0 sinh(n - i)/3o + -F„-i,o cosh(n - i)^o 
Enfi sinh(n + ^)/3o + cosh(n + i)/3o 


3 3 

En+ 1,0 sinh(n + -)/3o + -Fn+i.o cosh(n + -)/3o 


= 0 , (n > 1 ), 


(85) 


and 


sinh Po 


{2n - 1 ) 
+ cosh Po 
(n + l) 


Anfi sinh(n + i)/3o + cosh(n + i)/3o 

Cn-ifi sinh(n - i)/3o + -Dn- 1,0 cosh(n - i)/3o 


Cnfi sinh(n + i)/3o + D„,o cosh(n + i)^o 

3 3 

Cn+ifi sinh(n + -)/3o + -D„+i,o cosh(n + -)(5q 


(2n + 3) [ 


cosh/3o - 


2 n- 1 


„-(n-i)/3o _ ( ^ ^ ^ g-(”+§)/5o 

\ 2 tl “h 3 


( 86 ) 


(n > 0 ) 


Furthermore, Eqs. 46-48 and 51-52 can be used again in order to find the unknown coefficients, ^n, 0 : C'n ,0 5 ^n,o 
and Fn,o. 

By symmetry the force on the particle is in the ^-direction and there is no torque; explicitly 


CX) ^ 

= 2-\/^7r € Uz y; {Cnfi + Dnfl) F {n + -) {Anfi + Bnfi). 


(87) 


n=0 


E. Rotation of a non-translating sphere parallel to a solid wall 

This problem was first discussed by Dean and O’Neill [81].The boundary conditions on the wall are no-slip and 
no-penetration. 


K’1/3=0 ~ "^01/3=0 ~ "^2; 1/3=0 — O 5 


(88) 
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and the velocity at the particle surface is: 


Vr\p^ = {z-{5 + R))cos(l), (89) 

^^l/3o = “ + 

Vz\i 3 ^ =-r cos(f). (91) 

From the boundary conditions at the particle surface it can be shown that only the m = 1 term survives in the general 
solution of the Stokes equation and that 70 = 0 . 

By applying the boundary conditions and using appropriate relations among associated Legendre polynomials, we 
obtain recursion relation 72 along with 


(n - 1 ) 


2 ( 2 n - 1 ) 

{n + 1 ) (ti + 2 ) 


An- 1,1 sinh(n - i)^o + cosh(n - ^)f3o 


+ 


2(2n + 3) 
n 


3 3 

An+1,1 sinh(n + -)/3o + cosh(n + -)^o 


( 2 n - 1 ) 

+ cosh /3o 

(n + 1 ) 
“ (2n + 3) 

2 \/2e 


Gn-1,1 sinh(n - i)/3o + cosh(n - i)/3o 

Gn,i sinh(n + ^)/3o + cosh(n + i)/3o 

3 3 

Gn+i,! sinh(n + -)^o + Hn+i,m cosh(n + -)/3o 


e ("“'' 2 )^° + coth/3o 




^ sinh ydo 

for n > 0. Applying the boundary condition ( 91) into the expression (19), we have: 


sinh pQ 


An,i sinh(n + ]^)(5q + Bn,i cosh(n + i)/3o 


(^- 1 ) 
( 2 n - 1 ) 

cosh /do 

(n + 2 ) 
(2n + 3) 

V 2 e 


C„_i,i sinh(n - i)/3o + -Dn-i.i cosh(n - i)/3o 

C„,i sinh(n + i)/3o + -Dn,i cosh(n + i)/3o 

3 3 

Cn+ 1,1 sinh(n + -)^o + ^n+ 1,1 cosh(n + -)^o 


1 


2 n- 1 


3-(^-i)/5o 


+ 


277/ “ 1 “ 3 


^-(n+|)/3o 


, {n > 1 ) 


(92) 


(93) 


The no slip and no penetration condition at the wall yields Eq. 75-77. Two more relations obtained form mass 
conservation, Eq. 78-79. 

The unknown coefficients of velocity fields are found form simultaneous solution of above relations and the force 
and torque on the particle are [81] 


— —a/^tt e fi R Q.y (Gn,i+-f7n,i) ~ '^(^ + 1) (^n,l + ^ 71 , 1)5 


(94) 


n=0 


= 8tt ijR^ ny 


sinh^ /do 


00 


3 12^2 7, 

(tT/ (tT/ + 1) {2Cn,l “1“ coth -|- (277/ -|- 1 — coth /doGn,l)) 

(2 - j (n (n + 1 ) {2Dn,i + coth/3o-B„,i) + (2n + 1 - coth/3o-f7n,i)) 


(95) 


with all other components of force and torque vanishing. 
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